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' The spectrum of two-point functions in a holographic renormalization group flow from an ultraviolet 

(UV) to an infrared (IR) conformal fixed point is necessarily continuous. For a toy model, the spectral 
function does not only show the expected UV and IR behaviours, but other interesting features such as 
sharp peaks and oscillations in the UV. The spectral functions for the SU{3) x t/(l) flow in AdSi/CFTi 
Q^' and the SU{2) xU{l) flow in AdS^/CFTix are calculated numerically. They exhibit a simple cross-over 

^ ' behaviour and reproduce the conformal dimensions of the dual operators in the UV and IR conformal 

phases. 

(N ' 
> 
(N 

^ : 1 Introduction 



In the AdS/CFT correspondence, holographic renormalization group (RG) flows provide the dual gravita- 
tional description of strongly coupled gauge theories, the conformal invariance of which is broken either by 
an explicit operator deformation or by some non-zero operator vacuum expectation value. The calculation 
of correlation functions in such field theories via holographic renormalization [T] |2l O is an elegant and 



\-{ ' powerful method to obtain physically interesting quantities such as the mass spectra of glueballs |0] |5l and 

mesons (see IH for a review), scattering amplitudes EIH, as well as hydrodynamic transport coefficients 
of the high-temperature plasma phase |l9j[T0l[TTl. Holographic RG flows have been used as gravity duals of 
high-Jc superconductors |[l2lll3l[l4l[l5l and for constructing bottom-up AdS/QCD models lfT6l[T7l[T8l . Re- 
cently, progress has also been made in non-AdS/non-CFT cases, e.g., in the Klebanov-Strassler background 

nsiiioiEii. 

The present paper deals with holographic RG flows that interpolate between an ultraviolet (UV) and an 
infrared (IR) fixed point. In particular, we consider the gravity dual of the SU {2)xU{\) Leigh-Strassler flow 
ll22l|23]| and its three-dimensional cousin Il24l|25ll26l . which has been identified recently as a SU (3) x ?7(1) 
mass deformation of the Bagger-Lambert theory ||27l |28l |29l l30l . In both cases, the holographic RG flow 
involves two (active) scalar fields and is known only numerically. To our knowledge, correlation functions 
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in these backgrounds have not been studied so far. We intend to approach such calculations in the present 
paper by considering the two-point functions of the operators that are dual to the active scalars. In particular, 
we shall calculate the (eigenvalues of) the spectral function (matrix) and confirm that the IR physics of the 
conformal field theory is determined by the conformal operators living at the IR fixed point. 

The following simple argument shows that the spectrum of two-point functions in any RG flow to an IR 
fixed point with a gravity dual must be continuous. In AdS/CFT, the bulk dual of a mass state is obtained by 
imposing two conditions on the solution of the linearized bulk equations of motion, namely regularity in the 
bulk interior and absence of the dominant asymptotic mode (implying integrability). Typically, if there are 
regular and singular bulk solutions, these two conditions can be satisfied only for discrete mass values. If the 
deep bulk interior is described by an AdS geometry corresponding to the IR fixed point, the two deep bulk 
behaviours of a scalar field will be given by 

z-oo: 0,,,~z(^-i)/2e-^^S fcg-z^^-i'/^e'^^S (1.1) 

where d is the boundary dimension, L is the IR AdS radius, and k = {k^y/'^. For generic complex momentum 
we choose the square root such that Re^ > 0, so that the modes are regular or singular as indicated by 
the subscripts. However, for k^ = —m^ one has Re/c = 0, and both modes are wildly oscillating for large 
z. Hence, the regularity condition cannot be imposed, and the mass spectrum will be continuous with non- 
normalizable pseudo-states, as usual. 

Let us now outline the rest of the paper and summarize the results. In Sec. [2l we start by considering 
a toy model for a RG flow between an UV and an IR fixed point consisting of two patches of AdS space 
with different radii glued together. We choose the masses of the scalar fields such that the bulk solutions are 
simple. Moreover, the UV conformal dimension of the dual operator coincides with those in the AdS4/CFT3 
flow we study later. We verify that the spectral function exhibits the scaling behaviours in the UV (large 
mass) and the IR (small mass) that one expects for the operators living at the fixed points. In addition, we 
find interesting features in the spectral function such as sharp peaks and UV oscillations, which are are related 
to the discrete mass spectrum in the hard-wall model. Such features have been observed and interpreted as 
quasinormal modes in thermal spectral functions [31] . 

In Sec.[3l we present the technical tools that we need for calculating spectral functions in holographic RG 
flows. We start by reviewing the spectral function in many-particle systems, continue with a review of the 
equations governing the bulk dynamics in holographic RG flow backgrounds and then present the numerical 
strategy for the calculations. In Secs.HandU the SU{3>) x U{\) flow in AdS4/CFT3 and the SU{1) x U{1) 
flow in AdS5/CFT4 are considered, respectively. In both cases, we review the RG flow backgrounds and 
present our numerical calculations of the eigenvalues of the density function. Our results indicate that the 
deformed CFTs exhibit a simple cross-over behaviour from the UV to the IR regime. In the UV regime, the 
CFTs behave as the undeformed ones, while the IR behaviour is dominated by the operators that live at the 
IR fixed point. 
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2 Toy Model 



Let us consider the following simple toy model consisting of two patches of {d + 1) -dimensional AdS bulk 
space-time glued together at some fixed radial coordinate. We use an AdS metric of the form 

d^^ = dr^ + e^^*") riij dx' dx^ , (2.1) 

where 

A{r)=r/L+Ao, (2.2) 

L being the AdS length scale, and Aq an arbitrary constant. The equation of motion for a free massive scalar 
field can written in the form 



0=0. (2.3) 



The scalar's mass m is determined by the well known relation A = \/(P- /A + m^L'^, while the dimension 
of the dual operator is, in the case of regular boundary conditions, A = (i/2 + A. After introducing the 
dimensionless variables 

P = k^L'^e-^° , z = e.''''^ , (2.4) 
(12.31 1 takes a standard form and has the solutions 

0reg = Z''^ (kz) , 0sing = Z^^' Ia (fe) , (2-5) 

where K;^ Ia denote modified Bessel functions. Here and henceforth, k is defined as the square root of 
P with -n/2 < argk < n/2. 

In the following, let us consider a UV region, r > 0, with parameters Luv = 1 and Ao,uv = 0, in which 
lives a scalar field with Auv = 1/2- This choice coincides with the UV behaviour of the 3d RG flow that 
we study in Sec. H] For A = 1/2, the modified Bessel functions in (12.51 ) are elementary. Hence, the general 
solution of (12.31) becomes 



(rf-l)/2 



X 

scosh(fe) + - sinhte) 
k 



(2.6) 



where s and r are called the source and response coeffiencts, respectively. They determine the connected 
two-point function of the dual operator by 

in terms of which the spectral function p is defined as 

p = -2Im?^ , with P = -np- + /0+ . (2.8) 

In the IR region, r < 0, we set Lir = I, while the continuity of the metric at r = imposes Aqjr = 0. 
Let us choose Air = 1/2 + « for some positive integer n, so that the solutions ( 12.51 ) are again elementary. 
Imposing the regularity condition for r — oo (the deep interior) selects the solution 



3 



The overall numerical factor is irrelevant for what follows, as it drops out of ( I2.7I ). 

The source and response coefficients, s and r, are determined by appropriate matching conditions at 
r = 0, i.e., z = 1. From (12.31 ) and (12.41) follows that the matching conditions are 



z= 1 : ^ = 0iR , 
They give rise to 



d-l 



n 0IR 



1 / 



1 / r 
21^-^ 



n + j-\)\ 



T — ^ 



{n + j-\)\ 

f^jiin-mikiy 



so that one obtains, using ( 12.71 ). 



k + ', 



n 

J (7-1) !(«-;■) !(2W)^ 



j\{n-jy.{2kiy 



Curiously, the boundary dimension d has dropped out. Defining a by 



-2a 



— = — tanh a , 

A: 



and, similarly, y and j3 by 



(n+j-ry. 



U-\y.(n-jy.(2kiy 



l+n I 



(n+j-\y 



(I2l3l) yields 



^(k) 



1-y 



-2/3 



J = tanh j8 



where j is determined from ( 12.151 ) as 

\ {n-l+j)\ 



p,{n-\-j)\i\{lkl)j 



tanh/c + J 
1 + jtanh/c 



^ {n + i)\ 
p,{n-j)\i\{2kl)j 



-1 



K„- 1/2(^0 
K„+ 1/2(^0 



(2.10) 



(2.11) 
(2.12) 



(2.13) 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



In the following, let us extract the UV and IR behaviours of the spectral function and verify that they 
match with the spectral functions for operators of dimensions Auv = (<i + l)/2 and Air = {d + \)/2 + n, 
respectively. In the UV, i.e., for large \k\, y is approximately 3; 1 — n/{kl), such that ( 12.161 ) becomes 



UV : Sf(k) 



-k+-. 



2k 



(2.18) 
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For Rek > 0, only the first term on the right hand side remains in the UV limit, which is precisely the 
propagator of a A = 2 operator, as expected. In addition, the second term gives rise to oscillations in the 
spectral function (12.8b . which becomes 

UV: p(m2) w2Vw?+_sin(2\/n?) . (2.19) 

To obtain the IR (small m) behaviour of the spectral function, we start by inserting ( 12.161 ) into ( 12.81 ). 

p(„^) = -,M.+/)|[i^. = (2.20) 

Since, on the imaginary axis, k* = —k, only those terms of y that are even in k contribute to {y + y*). To find 
these, let us rewrite (12.171 ) as 

kl Pn-i{kl) 



2n-\ Pn{kl) ' 

where the polynomials P„(x) are defined by 



(2.21) 



" 1) ' 



It is straightforward to show that the P„ satisfy the recursion relation 



v-2 



Pn{x) = Pn-l{x) + -z r-TT rr/'n-l W (« > 2) , (2.23) 

[2n — \)[2n — 3) 



while the cases n = and « = 1 are simply Pq{x) = 1 and P\_{x) = \+ x, respectively. Hence, 

Pn-l{x) _ 1 

Pn{x) 1 + I Pn-2{x) ' 

^ ' ^+ (2n-l)(2n-3)P„-iW 



(2.24) 



which can be used recursively until one hits the term Po{x) /P\ (x). Therefore, expanding (12.241 ) for small x, 
the first odd term one encounters is the one that stems from Pq{x) /P\ (x) 1 — x. Thus, one has 

Pn-\(x) (-1)"(2«- l)x2"-l 

= even terms + -^^ + • • • , (2.25) 

P„[x) [(2?i-l)!!J^ 

so that (12.211 ) leads tcQ 

2(-l)"(W)2« 



= [(2.-i)np +-- ^'-'^^ 



Finally, inserting (12.261) into (12.201 ) yields 



, 2. 2/2«(m2)"+i/2 
IR: Pi^)= u... i\„l2 +••• • (2.27) 



[{In-r 



^This result can also be obtained, by similar reasoning, considering y in terms of the modified Bessel functions [see ( I2.17H and 
using the recursion relation Ky-i (x) — Ky+i (x) = — 2v/xKv(x). 
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The power behaviour m^"+' is expected for a dimension Ajr = {d + \)/2 + n operator. 

Another interesting limit is Z — > 0, which describes the set-up of a hard-wall cut-off. In this hmit, (12.211 ) 
and (12.221 ) imply j = 0, so that (12.161 ) becomes 



1 = 0: '^{\i) = -ktmhk. 



(2.28) 



In this case, the non-vanishing spectral function stems from the poles of tanhfc, which lie just on the imagi- 
nary axis. Setting k = im + e, one easily finds the discrete spectrum 



7=0 



(2.29) 



In the general case, the right hand side of (12.161 ) has an infinite number of poles on the left side of the k-plane, 
i.e., in the unphysical region Imk < oi The precise locations of the poles depend on n and /, but the limit 
/ ^ moves them onto the imaginary axis. For non-zero /, though, a number of poles tend to lie very close 
to the imaginary axis leading to a sharply peaked spectral function. Moreover, the UV oscillations in the 
spectral function ( 12.191 ) stem from the poles in the fc-plane. It is an interesting fact that these features have 
also been found in the context of thermal spectral functions at finite baryon density |[3T1 . The results are 
illustrated in Fig.[T] 



(a) 




10 15 20 25 30 35 40 
m 



(b) 




10 15 20 25 30 35 40 
m 



Figure 1: Spectral function p{m^) for parameters (a) « = 5, / = 4/5 and (h) n = 2,1 = \ /3. The sharp peaks 
are due to poles of the right hand side of (12.161 ) that lie in the unphysical region, but close to the imaginary 
axis k = im. The first peak of every plot is actually much taller, but the j-range has been cut off to show the 
other features. The UV behaviour ( 12.191 ) with oscillations is clearly visible. 



^Poles in the physical region ReA; > would give rise to poles of which in turn would be a signal of unstable states, or 

resonances. The absence of such poles implies the absence of resonances. 
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3 Spectral Functions in Holographic RG flows 



3.1 The Spectral Function 

The spectral representation of two-point functions is, of course, well known, but readers might remember 
only the simple case of a single operator. As the general, many-particle case will be needed later, we briefly 
review it in this section. 

Consider the two-point correlation function between operators and j. 



= y d''^^^Kx)^J(0))e"^'' , (3.1) 

and its spectral representation 



where Pij is the spectral function (matrix), often also called spectral density. For brevity, let us suppress the 
operator indices in what follows and regard ^^(k) and p{m^) as matrices. Considering complex k^, the two- 
point function may have poles and branch cuts along the negative real axis corresponding, respectively, to the 
discrete and continuous elements of the spectrum. The spectral function p is determined by the discontinuity 
of the correlation function across the branch cut, 

p(m2) = /(^+-^^_) , (3.3) 

where the subscripts + and — indicate the sign of the (infinitesimal) imaginary part of = —m^ + /0^. In 
addition, the correlation function (13.11 ) satisfies the analyticity property 

{k')]' = <;&[{k^y], (3.4) 



which implies 

Combining (13.51 ) with (13.31 ) yields 



'^l . (3.5) 



p(m2)=/(^+-^^j) , (3.6) 



which shows that p is a hermitian matrix with real eigenvalues. Eq. (13.61 ) reduces to the well known p = 
— 2Im^^+ for the case of a single operator. 

3.2 Bulk Dynamics 

Let us start by reviewing the equations governing the dynamics of the bulk fields |[32l[T9ll . which encode the 
information about two-point functions in holographic renormalization group flows. Then, we will introduce 
a change of variables that further enhances the system eliminating the need to calculate the warp factor. 
The systems we consider are of fake SUGRA type with actions of the form 

1 1 



+ Sh, (3.7) 
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where the potential V{^) is given in terms of a superpotential W{^) by 



V{^)=^-G"'WaWt-^ 



(3.8) 



We will not specify the boundary terms Sb in (I3.7K as they do not affect the bulk dynamics, although they 
are important for holographic renormalization. Our notation agrees with |fT9]| . In particular, field indices 
are covariantly lowered and raised with the sigma-model metric Gat and its inverse, G"^, respectively; Wa = 
daW = dW{<p)/d<p"; and covariant derivatives with respect to the fields are indicated by Da or by a "|" 
preceding the index, as in Wa\t = DtWa = dtWa — ^^/j^c, being the Christoffel symbol for the metric 

Gab- 

Holographic renormalization group flows are described by domain wall backgrounds of the form 



d^^ = dr^ + e^'-'-'^ 7]ij dx' dxj , 0" = 0" (r) , 
which satisfy the BPS equations 



drA 



d-l 



(3.9) 



(3.10) 



Linearized fluctuations around the domain wall background are best described in a gauge invariant fash- 
ion, in which the relevant, independent fields are the scalar fluctuations a" and the traceless transversal metric 
fluctuations, e'j. They satisfy the (linearized) equations of motion 



2d 



and 



D,- + M-- — jWj (D,-M) +e"^'*n 
2d 







d-l 



W 5, + e-^n 







respectively. In (13.111 ). we have omitted the field indices, M denotes the matrix 

WWb 



M\ = W"\b 
and Dr is the background covariant derivative 



W 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



For more details, we refer the reader to the original papers ir71 [32l[T9]| . 

Let us now introduce a change of variables, which facilitates the treatment of the bulk dynamics by 
eliminating the warp function A{r). Introducing a new radial variable a by 



:A(r) 



(3.15) 



where we implicitly assume that A (r) is a monotonous function, or equivalently, that W does not change sign 
along the RG flow, the BPS equations (13.101) yield simply 

d-\W"{^) 



Bar 



2 W{^) 



(3.16) 
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and the field equations ( I3.11I ) and ( I3.12I) become 



{Da+M + d) {Da -M) + e-2^ ^^-^U 



4W2 



= 



and 



e;. = o 



respectively. In (13.171 ). we have 



and 



d-\(W"\, W'^Wt\ d-\ (w- 



m 

As an aside, we observe that da^^ is a solution to the linearized field equation (I3.17I) . 



(3.17) 



(3.18) 



(3.19) 



(3.20) 



3.3 Near Boundary Asymptotics and Pure AdS Case 

To extract the CFT data from the bulk dynamics, we make use of holographic renormalization 111121 [3l- The 
near boundary behaviour (large a) of the solutions of (13.171 ) is determined by the UV fixed point, at which 
the matrix M is chosen diagonal with eigenvalues A,- — d/2. This follows from the expansion of W around 
the fixed point, which for a single scalar ^ reads lEl 



d-l 1 
' 2L ~4L 



(J-2A)02 + -- - . 



(3.21) 



The (UV) conformal dimensions of the dual operators are related to the eigenvalues of M by A,- = J/2 + 
Notice that the A,- can be negative. 
A generic solution of ( 13.171 ) can be decomposed as 



a''(a)=s;af(a)+r,a?(a) 



(3.22) 



where a, and a, denote the dominant and sub-dominant series solutions, which are the duals of the CFT 
operators ^, of dimensions A,, and 5, and t; are the source and response coefficients, respectively!^ The 
leading terms of o; and a; agree with the respective counterparts in pure AdsB 



—1— e-(^-^-> 5f + • • • , a^'(a) = e"^'^ 5f + • • • . 
2Ai-d 



(3.23) 



The somewhat unconventional normalization of the dominant solutions eliminates the proportionahty factor 
2 1 A, I in the standard relation between the exact one-point function and the response coefficients Il33ll34l . The 
sub-leading terms in ( 13.231 ) depend on the background and must be considered case by case. 



^Note the change of notation with respect to (2), where the source and response coefficients were denoted by'and", respectively. 
'^In the special case A; = d/2, one has o = pe^''''/^ H . 
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To calculate the Green's function, one imposes a regularity condition on a. If a has components, then 
there are rig independent regular solutions. Hence, s and r can be considered as rig x ns matrices, the first 
index labelling the regular solutions and the second index the asymptotic (sub-)dominant ones. Then, the 
Green's function (matrix) has the elegant form 

W = 5-^-x, (3.24) 

where the matrix multiplication sums over the index labelling the regular solutions. 

For a scalar field in pure AdS with generic mass (A non-integer), one has the standard result f35\ 



n-\M) 



2\M 



and the spectral function (13.61 ) i 



P('«) = TF7^^WfTy'" • (3-26) 



[r(i + |A|)]n2 



3.4 Numerical Strategy 



The numerical strategy we employ for calculating two-point correlation functions and their spectral functions 
in holographic renormalization group flows between two conformal fixed points closely follows what one 
would do for an analytic calculation. The starting point is the equation of motion (13.171 ). which we rewrite 
as a system of first-order ordinary differential equations (ODEs), 




M I 




(3.27) 



where some auxiliary components b have been introduced, and I is a unit matrix. In this paper, we deal 
with canonical scalar fields, so that the derivative is simply D(j = da- As we must consider complex = 
— m^ + ZO"'", with some tiny imaginary partj^l we must also treat the real and complex components of the 
scalars as independent, so that, for Us (complex) scalar fields, ( 13.271 ) is a system consisting of ODEs. 
Hence, decomposing each scalar field as 

yim I 

the complex momentum k becomes a 2 x 2 matrix, 

— \mk\ 



(3.28) 



(3.29) 

Rek ) 

We are interested in the solutions of (13.271 ) that are regular in the bulk interior, i.e., which behave as 0reg 
in (11.11 ) to leading order. However, using the exponential damping to impose this behaviour is not practical. 




^In contrast to l l3.25t . formula l l3.26t is valid also for integer A 
^As we mentioned in the introduction, an IR AdS 
distinguishing between regular and singular solutions 



^As we mentioned in the introduction, an IR AdS region in the bulk interior implies that Re A: > is a necessary condition for 



10 



because Rek is tiny, and the bulk solution typically undergoes wild fluctuations in the region dominated by 
the factor before reaching the physically interesting domain wall. Thus, it is better to define the 

fields 

which satisfy the field equation 

We numerically integrate (13.311) for a < Go and then continue with (13.271 ) for a > Gq, with Gq chosen in 
vicinity of the domain wall. The initial condition that we impose in the IR can be determined from the 
leading and next- to leading behaviour of a', which is 

a^-oo: a' = e-^^(l + 7e^ + ---)o[) , (3.32) 

where Oq is a constant field vector, and the matrix 7 is determined from (13.311 ) as 



k* 



2LiR|^r 



1 , . ' 

-+DaM+{M + d)M 



(3.33) 



4 

There are Ifig independent initial conditions, which give rise to 2ns independent regular solutions. The 
field redefinitions (13.301 ) have been chosen such that the vector (o'b') matches smoothly with (ab) at a = 
Gq, but omitting the term containing Gq would just amount to a multiplication by an irrelevant (complex) 
normalization factor. 

In the UV region, the source and response coefficients, s and r, respectively, can be extracted from the 
numerical solutions by using the decomposition (13.221) into dominant and sub-dominant terms. Then, the 
Green's function is obtained from (13.241 ) as a 2ns x 2«.v matrix consisting of ns x blocks of 2 x 2 matrices 
of the form 

/ Re% Im%\ 

Finally, (13.61 ) yields the spectral function (matrix) p, the eigenvalues of which are then calculated. 

To test the numerical strategy described above, we have numerically calculated the spectral density for 
the toy model studied in Sec.|2l In particular, we have considered the parameters d = 3, Ljr = 0.8 and n = 2, 
i.e., the scalar in the IR region is dual to an operator of dimension Air = 4. The integration range was chosen 
as — 4 < r < 8, and the IR and UV regions are glued together at r = 0. 

The exact solution for the spectral function is given by (12.81 ) and (12.161) . The results are depicted in Fig.[2j 
and the agreement between the exact and the numerical solutions is evident. 



4 SU{3) X U{1) RG flow in J = 3 
4.1 Fixed Points and Background Solution 

We wish to study the case of the SU{3)xU{l) invariant RG flow in AdS4/CFT3. The ^ = 2 domain wall 
solution in ^ = 8, 4-dimensional gauged SUGRA was found in Il24ll25]| and Ufted to M-theory in ll26l . 
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Figure 2: Spectral density of the toy model with parameters d = 3, Ljr = 0.8, n = 2. The numerical solution 
(points) agrees with the exact solution (line). 



and the superpotential 



The fake SUGRA system describing the flow contains two scalar fields, p and x, with the sigma model 
metric 

Gpp = \2p-\ G^^ = \, Gxp=0 (4.1) 

W = - ^ [cosh(2;t) + 1] + [cosh(2;t) - 3] . (4.2) 
It will be helpful to introduce also a canonical basis {cc,x) of scalars by defining 

p = e5^ . (4.3) 

The background equations (13.161) read 



16W 



-l^ (cosh(2;t;) + 1) + p^(cosh(2;t;) - 3) 



5.Z = -^(p^-^)sinh(2^) 



(4.4) 



The fixed points we are interested in are the (attractive) UV fixed point at p = 1 (a = 0), ;^ = and the IR 
fixed point at p^ = 3, cosh(2;^) = 2. (Notice that p is assumed positive, and amongst the two equivalent 
IR fixed points we pick the one with positive X-) The solution of (14.41 ). which interpolates between the 
fixed points, can be found numerically and is illustrated in Fig.|3] For numerical stability, one sets the initial 
conditions for the field values very close to the IR fixed point, and the dynamics of (14.41 ) lets the solution flow 
to the UV fixed point with increasing a. This numerical solution will serve as the domain wall background 
in the numerical treatment of the fluctuations, which follows. 



^We have adopted the expressions of 1241 1251 126| to the conventions of Sec. l3.2l and set the length scale of the asymptotic AdS 
region in the UV to unity. 
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Figure 3: (a) Contour plot of W(p,;^) with the solution of (14.41 ) connecting the UV and IR fixed points, (b) 
Plot of W vs. a along the flow. (c,d) Plots of the scalar fields vs. a. 

Before continuing, we would like to make some comments on the integration constants for ( 14.41 ) Fl For 
large values of a, i.e., in the UV region, the scalar fields approach their respective UV fixed point values as 



large a 



a{o) 



^ —(7 

ae 



(4.5) 



with two coefficients a and %■ However, only their ratio, or equivalently, an angle in field space, determines 
the direction of the flow. For the flow starting at the IR fixed point, the numerical solution yields 



a 

- f« 0.6529 . 
X 



(4.6) 



A possible common factor can be absorbed by a shift of the radial variable o ^ o' = o + 8o . 

Let us now consider the dimensions of the dual operators at the fixed points. This was done in ll24l . but 
we shall rederive them here for completeness. The dimensions of the operators dual to the scalars a and x 
are determined by expanding W about a fixed point to quadratic order. At the UV fixed point, one finds 



(4.7) 



similar discussion for tiie flow in AdS5/CFT4 can be found in 1221 and in Sec. 15. II 
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Comparing this to the generic formula ( I3.21I ). we find that the scalars a and % correspond to a doublet of 
relevant operators of dimension A = 2. 
For the IR fixed point we introduce 

a = aiR + a, X = Xir + X, (4.8) 
where aiR and ;^ir are the respective field values at the fixed point. Expanding W to quadratic order yields 



1 



^IR 



1 



W = —— l + -a^-lax + 



with LiR = 2 • 1>^^I^. Rotating the fields by 



where 



cos 9 



a' 

X 



in. 



cos (p —sm(p 
.sinq) cosq) 



1 



1/2 



a 



1/2 



one brings the quadratic terms in ( 14.91 ) into diagonal form and obtains 
Hence, the dual operators and have dimensions 

1 / , N 

A 



1 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



(4.13) 



From the monotonicity relation drW = WaW" > 0, which stems from ( 13.101 ). follows that the flow approaches 
the IR fixed point along the x' direction, and the irrelevant operator i^^/ controls the RG flow in the field 
theory 1,22 J . 



4.2 Spectral Functions 

Here, we present our numerical results for the spectral functions of the two-point correlators for the dou- 
blet of operators 0'a and i^^,, which have UV conformal dimensions A^ = A;^ = 2. As shown above, the 
renormalization group flow lifts the degeneracy of the dimensions and ends at the IR fixed point with two 
operators and ^^z, the dimensions of which are given by (14.131 ). From the eigenvalues of the spectral 
function matrix, we shall be able to identify clearly the cross-over from the UV to the IR. Furthermore, we 
shall look for evidence of UV oscillations in the spectral functions, as in the toy model, but find no evidence 
of them. 

The system of field equations (13.271 ). which we have to solve numerically, contains in o the linearized 
fluctuations of the scalars a and X- As the kinetic term for these scalars is canonical, we have simply 
Da = da- Furthermore, the matrix M takes the form 

M= — ( T[cosh(2;^)-3][cosh(2z) + l] ^p4sinh(2x) 

W2 1^ f p4sinh(2z) j^(p«-3)[3(p« + l)cosh(2;^)-p« + 3] 
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where the superpotential W is given by p = e«/P%/3)^ jhe fields are evaluated on the (a-dependent) 
background. As explained in Sec. 13.41 we use (13.311) in the IR region for numerical stability choosing Gq = 0. 

The results for the eigenvalues of the spectral function matrix are shown in Fig. |4] and exhibit the fol- 
lowing features. For large m, both eigenvalues show the expected UV behaviour ( 13.261) with A = 1/2, 
i.e., p = 2m, without any sign of oscillations. (We have checked also for higher values of m than those 
shown in the figure.) For small m, a fit of the log of the larger eigenvalue vs. Inm yields the formula 
Inpi 2. 12 Inm + 1.34, and the coefficient of In/n agrees nicely with the expected value for the operator 



2Aa/-3 = \^-2«2.12 



(4.15) 



A similar fit for the smaller eigenvalue in the IR region yields lnp2 « 6. 17 In m + 1.21. The slight disagree- 
ment of the slope from the expected value, which is 6. 12, is due to the larger numerical error for 

The results for the spectral functions indicate that the deformed CFT exhibits a simple cross-over be- 
haviour under the RG flow, with a cross-over point m = A^. For m > Ac, its dynamics is governed by the 
undeformed ^ = 8 CFT, whereas for m < A^. the IR conformal phase takes over. The numerical value of 
Ac is subject to shifts of the radial variable a ^ a' = a + So, which just amounts to choosing a momentum 
scale. 



100 



(a) 




100 



(b) 




Figure 4: Log-log plots of the spectral function eigenvalues pi (a) and p2 (b) vs. m. The straight lines 
correspond to the expected IR and UV behaviours. The UV lines are given by p = 2m, while the IR hnes are 
(a) pi = 3.83m2^2 (b) p2 = 2.Sm^-^^. 



5 SU{2) X RG flow in J = 4 
5.1 Fixed Points and Background Solution 

This RG flow was discussed extensivley in ll22l . Here, we will give a brief review for completeness. 



'The smaller eigenvalue P2, which is several orders of magnitude smaller than pi in the region of interest, is calculated as the 
difference of two numbers of the order of pi . 
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The fake SUGRA system describing tiie RG flow contains two scalar fields, p and with the sigma 
model metric 



Gpp=6p-\ G^^ = l, G^p=0 (5.1) 



W = -2^[cosh(2;t) + l] + ip^[cosh(2;t)-3] . (5.2) 



and the superpotential 



A canonical basis of scalars j8 , ;^ is achieved by introducing 

p =e^ . (5.3) 

The background equations (13.161) read 



^(cosh(2;t:) + 1) + p4(cosh(2;^) - 3) 
P 



5.Z = -^(/-^)sinh(2;,) 



(5.4) 



The fixed points we are interested in are the (attractive) UV fixed point at p = 1 (j8 = 0), ;^ = and the 
IR fixed point at p^ = 2, cosh(2;^) = 5/3. The numerical solution of (I5.4I ). which interpolates between the 
fixed points, is illustrated in Fig. [5] This numerical solution will serve as the domain wall background in the 
numerical treatment of the fluctuations, which follows. 

The background behaviour of the fields near the UV fixed point is 

largea: zW^^e"'^, i8(r) p. jSe-^^+y^f ae-^^^ , (5.5) 

with two coefficients j8 and %■ The combination of them, which is invariant under a shift of the radial 
variable, a ^ a' = a + 8o, is _ 

Z' ' V 3 
which agrees with the result of 1221 . 

Let us rederive the dimensions of the dual operators at the fixed points. At the UV fixed point, one finds 



^ + \/:7lnZ~ -1-469(6) , (5.6) 



W = -{^- + ^' + \x^^+- . (5.7) 

Comparing this to the generic formula (13.21b . we find that the scalars j3 and % correspond to relevant opera- 
tors of dimensions ^^=2 and A_;j. = 3, respectively. 
For the IR fixed point we introduce 

i3=j8iR + j8, Z = ZiR + Z, (5.8) 
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In the notation of (22|, P = (pj, X = 'P\- 
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Figure 5: (a) Contour plot of W{p,x) with the solution of (15.41 ) connecting the UV and IR fixed points, (b) 
Plot of W vs. a along the flow. (c,d) Plots of the scalar fields vs. a. 

where /3ir and Xir are the respective field values at the fixed point. Expanding W to quadratic order yields 

1 / 3 



W 



with LiR = 3 • 2"^/^ Rotating the fields by 

I P'\ /cos(p — sin^A /j8\ 



where 



cos<p 



1 f 1 



.1/2 



sm(p 



1 / 1 



1/2 



one brings the quadratic terms in (15.91 ) into diagonal form and obtains 

?ip, = l-V7 , A^/ = 1 + . 

Hence, the dual operators and have dimensions 

Ap, = l + Vl, Ajf/ = 3 + \/7. 
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(5.9) 



(5.10) 



(5.11) 



(5.12) 



(5.13) 



The background approaches the IR fixed point along the %' direction, and the irrelevant operator controls 
the RG flow in the field theory. 

5.2 Spectral Functions 

Here, we present our numerical results for the spectral functions of the two-point correlators for the operators 
and which have UV conformal dimensions A^g = 2 and A;^ = 3, respectively. As shown above, the 

renormalization group flow ends at the IR fixed point with two operators 6^, and G^i, the dimensions of 

which are given by (15.131 ). As in the AdS4/CFT3 case, the spectra show a cross-over behaviour along the 

flow and no signs of oscillations. 

The system of field equations ( 13.271 ). which we have to solve numerically, contains in o the linearized 

fluctuations of the scalars j8 and %. As the kinetic term for these scalars is canonical, we have simply 

D(j = ^cj- Furthermore, the matrix M takes the form 

J_ /¥['^osh(2;^)-3][cosh(2z) + l] ^flp^ ix^^x) 

^p2sinh(2;t;) (p^ - 2) [(3p^ + 2)cosh(2z) - p^ + 2] 

where the superpotential W is given by (Eal), p = e^^/^^ and the fields are evaluated on the (a-dependent) 
background. 

The results for the eigenvalues of the spectral function matrix are shown in Fig. |6] and exhibit a clear 
cross-over behaviour from the UV to the IR. For large m, the spectral functions exhibit the expected UV 
behaviour (13.261) for operators of dimensions A^ = 2 and ts.^ = 3, 

Pi=pp?s7r, p2 = p;^ = |m^ (5.15) 

respectively. Again, there is no sign of oscillations. 

For small m, a fit of Inp vs. Inni yields the relations pi w 3.51nm + 4.3 and p2 ~ 7.61nm + 3.5, while 
the values expected for the slopes are, respectively, 

2(\/7-l) «3.29 , 2(^7 + 1) ?s 7.29 , (5.16) 

as one obtains from the IR dimensions ( 15.131 ). We suspect that the slight disagreement stems from numerical 
issues, but this question deserves further investigation. We just note that, with a numerical accuracy of the 
integration routine of about 10^^, it is impossible to calculate the eigenvalues for lower m. At m = 0.08, the 
eigenvalues differ already by five orders of magnitude. 
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Figure 6: (a) and (b) log-log plots of the spectral function eigenvalues pi and p2 vs. m, respectively. The 
straight lines correspond to the expected IR and UV behaviours. The UV lines follow from (13.26b with (a) 
A = and (b) A = 1. The IR hnes are (a) pi = 40m3 ^9 ^ l5m^-29. 
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